The impact of nonlocal response on metallo-dielectric multilayers and optical patch 

antennas 
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We analyze the impact of nonlocality on the waveguide modes of metallo-dielectric multilayers 
and optical patch antennas, the latter formed from metal strips closely spaced above a metallic 
plane. We model both the nonlocal effects associated with the conduction electrons of the metal, 
as well as the previously overlooked response of bound electrons. We show that the fundamental 
mode of a metal-dielectric-metal waveguide, sometimes called the gap-plasmon, is very sensitive to 
nonlocality when the insulating, dielectric layers are thinner than 5 nm. We suggest that optical 
patch antennas, which can easily be fabricated with controlled dielectric spacer layers and can be 
interrogated using far-held scattering, can enable the measurement of nonlocality in metals with 
good accuracy. 



I. INTRODUCTION 

With the emergence of new analytical, numerical and 
nanofabrication tools, the pursuit of plasmonic systems 
for a variety of nanophotonic applications has expanded 
rapidly in recent yearsi~— . Plasmonic media here can be 
defined as conducting surfaces and nanostructures, whose 
optical scattering is largely dominated by the response of 
the conduction electrons. Plasmonic behavior is typically 
associated with excitation wavelengths at which the in- 
ertial inductance of the charge carriers plays a critical 
role in the collective response -. In the design of plas- 
monic media, the dynamics of the conduction electrons 
can often be well approximated by assuming a Drudc- 
like model for the permittivity, which has the frequency 
dispersive form 



e = 1 - 



uj- 



iTuj 



(1) 



assuming a time dependence of e~ lut . The plasma fre- 
quency, uj p , proportional to the square root of the carrier 
density, typically lies within the ultraviolet portion of the 
spectrum for many metals. Thus, for frequencies just be- 
low the plasma frequency, the electric permittivity can 
be characterized as a lossy dielectric, for which the real 
part of the permittivity is moderately negative. At wave- 
lengths where the real part of the permittivity is nega- 
tive, surface plasmon modes can be supported, which are 
collective oscillations of the coupled electromagnetic field 
and conduction electrons. Surface plasmons can serve to 
transport energy along metal surfaces in a manner similar 
to dielectric waveguides, but are also playing an increas- 
ingly important role in the field of metamaterials, where 
metallic nanostructures are often used as elements that 
provide strong and customizable scattering. Surface plas- 
mons represent the underlying mechanism behind perfect 
lenses^, hyperlenses^, spasers^— and many other pro- 
posed metamaterial-related devices. 

The simplicity of the Drude model of electron response, 



([T]), has has enabled the rapid modeling of plasmonic and 
metamaterial structures; The salient features associated 
with most plasmonic structures presented to date can 
usually be computed with sufficient accuracy -sometimes 
even analytically- assuming the Drude formula. Par- 
ticularly when the underlying physics is the main focus 
rather than detailed performance characteristics, Eq. ([1]) 
frequently provides an adequate description of the plas- 
monic response. It should be noted that despite the 
relatively simple form of Eq. (JXJ) , the numerical simu- 
lation of plasmonic systems remains a non-trivial task 
because the surface plasmon spatial variation is not lim- 
ited by the wavelength of light; rather, the surface plas- 
mon can confine light to nanometer sized regions, mak- 
ing plasmonic structures an inherently multiscale model- 
ing problem^ 2 -. Thus, the frequency dispersion and the 
negative permittivity associated with the Drude model 
contain non-trivial physics, and have been successfully 
applied to a wide range of plasmonic and metamaterial 
configurations. Naturally, the actual electronic response 
of a metal or highly doped semiconductor is much more 
complicated than that suggested by Eq. (1). Plasmonic 
structures are now reliably fabricated at the nanome- 
ter and sub-nanometer scales, where new optical prop- 
erties arise that cannot be accounted for solely by the 
Drude mode l 13 i 14 . Since these sub-nanometer features 
are likely to be crucial for optimizing field localization 
and enhancement^ - — , a more detailed description of the 
properties of plasmonic devices is demanded. Effects that 
would be secondary or of no consequence to the overall 
function of prior plasmonic devices, may introduce ma- 
jor constraints on the detailed performance and ultimate 
competitiveness of optimized plasmonic structures with 
subnanometer features. For these reasons, it is relevant 
to consider a more advanced physical model of the carrier 
response in conductors. 

The Drude model of a conductor assumes only the par- 
ticipation of conduction electrons (no bound charges), 
and further assumes a straightforward force-response re- 
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lationship between the applied electric field and respond- 
ing current density. An intrinsic feature of this model 
is that the responding current density at a given point 
within the material is proportional to the electric field 
at that point; that is, the Drude model assumes local- 
ity. Even when the exact, measured values of the bulk 
permittivity are used, there is an implicit assumption of 
locality since the permittivity is only a function of fre- 
quency rather than of both frequency and wave vector. 
To capture the additional physics associated with elec- 
tronic response, it is necessary to consider a more de- 
tailed model of the force-response relationship between 
the field and current density. 

More accurate descriptions of the free electron gas 
have been proposed in the past, including a description 
based on a hydrodynamical model for the conduction 
electrons 1 ^—, and a microscopic description initiated by 
Feibelman 26-2 -. The latter has been improved over the 
years 29 and has been recently used to include the effects 
of nonlocality on metallic slabs2£ and slot waveguides^. 
The hydrodynamical approach clearly suffers from an un- 
certainty about which additional boundary conditions 
should be used, but allows for more transparent phys- 
ical interpretations 2 ^. Moreover, the hydrodynamical 
model can be reasonably implemented in numerical cal- 
culations, and also is useful for finding closed-form, ana- 
lytical results. For example, the hydrodynamical model 
has been used in conjunction with transformation optics 
techniques to find analytical expressions for nanostruc- 
tures that illustrate the impact of nonlocal respons e 14 ' 15 . 
Recent experiments have shown that the hydrodynami- 
cal model is able to describe very accurately the plasmon 
resonance shift exhibited by spherical nanoparticles inter- 
acting with a metallic filmic. While the hydrodynamical 
model is clearly not the most sophisticated approach to 
describe the free electron gas, it can obviously capture the 
physics of nonlocality and it seems it can be made quite 
accurate through a correct choice of the free parameters 
it contains for situations of interest in plasmonics. 

In the present work, we first try to describe the re- 
sponse of the bound electrons as a polarizable medium, as 
has been shown to be accurate for Feibelman's method 29 , 
but in the framework of the hydrodynamical model. We 
find that this description greatly simplifies the discussion 
with respect to the additional boundary conditions. We 
explore the consequences of the model on the reflection 
of a wave by a metallic surface, on the surface plasmon 
and finally on the propagation of a guided wave along a 
thin metallic waveguide, as Wang and Kempa have shown 
that nonlocal effects could be expected^! for such a struc- 
ture. Using an analytical dispersion relation, we show 
that the nonlocal effects are enhanced in the slow light 
regime, when the waveguide is a few nanometers thick. 
Finally we study the large impact of nonlocality on opti- 
cal patch nanoantennas^ 2 - - — where the gap beneath the 
patch behaves as a cavity, making these structures ex- 
tremely sensitive^. The optical patch geometry paves 
the way for future experiments in which the effects asso- 



ciated with nonlocality will have easily measurable effects 
at wavelengths in the visible. 

II. NONLOCAL RESPONSE OF METALS 

While our analysis is not specific to metals, we use the 
term metal throughout while keeping in mind the analy- 
sis can be applied to highly doped semiconductors^! and 
potentially other conducting systems^. The polariza- 
tion of a metal, and hence its dielectric function, gen- 
erally contains contributions from both bound and free 
conduction electrons. Because we need to apply differ- 
ent physical response models to the free and bound elec- 
trons, it is essential to first distinguish their relative con- 
tributions. The experimental permittivity curves can be 
fit^S with a Drude term ([TJ that models the free electron 
contribution, to which is added a sum over the Brendel- 
Bormann 4 ^ oscillator terms that models the susceptibil- 
ity arising from the bound electron contributions. Fig- 
ure [1] shows the permittivity of gold obtained through 
the model as well as the fitted Drude permittivity, cor- 
responding to 1 + x/, where \f is the the susceptibility 
of the free electrons. The difference (not shown) between 
the modeled permittivity and the fitted Drude term cor- 
responds to the contribution of the bound electrons, \b- 

We assume here that the nonlocal response of the metal 
is largely dominated by the nonlocality induced by free 
electrons, so that we can treat the bound electron contri- 
bution as purely local, as some authors do 14 ' 41 ' 42 . Bound 
electrons too can be expected to present a nonlocal re- 
sponse, similar to what occurs in dielectric o 43 ' 44 . Our 
assumption is equivalent to assuming that the interac- 
tions between electrons in a free electron gas (through a 
quantum pressure and Coulomb repulsion) are more in- 
tense than essentially dipole-dipole interaction between 
bound electrons. 

Under an applied electric field, the medium will un- 
dergo a polarization with contributions from both bound 
and free electrons. The total polarization vector can thus 
be written 

P = P b + P / (2) 

where P^ = eoXfcE, \b being the susceptibility of the 
bound electrons, with the currents in the free electron 
gas related to the polarization in the usual manner: 

P/ = J- (3) 

By incorporating all responding currents and charges 
into the polarization, we can treat the metal as a dielec- 
tric, such that the electric flux density, D = eoE + Ph + 
P/, satisfies V • D = 0. Taking the divergence of D and 
writing P^ in terms of the electric field, we obtain 

V-P/ = -e (l + Xfc)V-E. (4) 

where we have explicitly assumed that the bound elec- 
tron susceptibility is local and can be taken outside the 
divergence operator. 
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FIG. 1. Real part of the relative permittivity of gold in the 
visible (solid line) and Drude permittivity according to— (dot- 
ted line). The difference is the contribution of bound elec- 
trons. Inset : imaginary part of the relative permittivity, 
same x scale, same lines. 



The free electron current density can be related to the 
applied electric field using the hydrodynamical model. 
Using Eq. 2, a linearized equation relating P/ to the 
electric field is given by 2 ^ 



-/3 2 V(V-P/)+P/+7P/ =eowj|E 



(5) 



where 7 is the damping factor, due to collisions of the 
electron gas with the ion grid, uj p is the plasma frequency 
of the metal, and (3 is the phenomenological nonlocal pa- 
rameter, proportional to the Fermi velocity vp . Usually 

the value of (3 = \J\vf has been considered in the liter- 
ature. However, a slightly more realistic hydrodynamic 
model should take into account other sources of nonlocal- 
ity, such as the Bohm potential, which can be shown to 
be of the same order of the Fermi pressure^. Though it 
is beyond the scope of this paper to introduce a more so- 
phisticated model, it makes sense from a phenomenolog- 
ical approach to consider a more empirical value for the 
parameter /3. Recently it has been shown for plasmonic 
systems of film-coupled gold nanoparticles that the value 

(3 — \J ~ 1.27 x 10 6 m/s gives a very good agree- 
ment with experimental datai^. In this work we will then 
assume this former value for both gold and silver. 

Assuming a harmonic solution of the form e~ lLJt , and 
using equation 21 the polarization Pf can finally be writ- 
ten 



-e - 



where the term 



e-(i + x 6 )^v(v-e; 



Xf 



(6) 



(7) 



can be identified as the local susceptibility associated 
with free electrons, corresponding to the Drude model. 



We have written the polarization terms in such a man- 
ner that the free and bound electron contributions can 
be distinguished. In determining the various parameters 
in these equations for the calculations that follow, we use 
the model provided by-22 and shown Fig. 1. 



III. TRANSVERSE AND LONGITUDINAL 
MODES IN METALS 

In a metal, taking the above description of nonlocality 
into account, Maxwell's equations can now be written 

V x E = iuj f i a H (8) 
VxH=-^(e (l + Xfc)E + P/) (9) 
= -iuse e (E - aV (V ■ E)) (10) 

where e is the local relative permittivity of the metal 



and 



1 + Xb + Xf 



Xf(l + Xb) P 2 



(11) 

(12) 
(13) 



l+Xb 



As shown rigorously in the appendix, there are two 
different solutions to these equations corresponding to 
two different kinds of waves. The first solution satis- 
fies V • E = 0, so that it corresponds to the standard 
solution to Maxwell's equations when the nonlocality is 
overlooked. Equations (j5]) and (|10p become 



V x E = zcj^oH 
VxH = — iwe eE, 



(14) 
(15) 



Finally all the fields satisfy Helmholtz's equation 

V 2 H + efcgH = 0, (16) 

2 

where fco = 72- ■ Since the divergence of the electric field 
is zero, the electric field is orthogonal to the wavevector 
when the wave is propagative, which means it is trans- 
verse. The dispersion relation for these transverse waves 
is thus 



(17) 



The second kind of solution is curl free, which means 
it satisfies V x E = and there is no accompanying mag- 
netic field. These waves are called longitudinal because 
when they are propagative, the electric field is parallel 
to the wavevector. They correspond to bulk plasmons: 
oscillations of the free electron gas due to the pressure 
term. Since the divergence of the electric field is not 
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identically zero, there exists a charge density inside the 
metal given by 

p = e V-E. (18) 

Equation (fT0|) . then yields the wave equation for the bulk 
plasmons 

V (V • E) - — E = V 2 E - — E = 0. (19) 
a a 

and the corresponding dispersion relation is 

k 2 = -- = i ^ 2 --^- + z 7W V (20) 
a P 2 \ 1 + Xfc j 

An alternative way to write this dispersion relation is 

n = i + x6- 2 , . p R2 , 2 = o, (2i) 

which is the way previous works have taken \b into 
account 1 - 4 " through a so-called longitudinal permittivity. 
But the equation governing the polarization (equa- 
tion (|5|)) cannot be deduced from the longitudinal per- 
mittivity using a simple Fourier transfor m 41 ! 42 , as has 
been previously pointed oul4£. 
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FIG. 2. Dispersion relation for transverse waves (thick lines) 
and bulk plasmons (thin lines) when (a) the nonlocality is ab- 
sent (fi = 0) and (b) the nonlocality is important. The dashed 
lines show the imaginary parts of the propagation constants 
for the transverse and longitudinal waves, below the plasma 
frequency oj p . When nonlocality is present, for a given ui two 
waves must be taken into account. (color online) 

The dispersion relations, Eqs. (fl~7|) and ([20]) . are plot- 
ted in Fig. [5] for two cases of the nonlocal parameter (3, 

for the simplified case where e = 1 — — §. For small j3, the 
longitudinal mode disperses very little, and can be gen- 
erally ignored in wave propagation problems. When (3 
is nonzero, however, the longitudinal mode acquires dis- 
persion, and is generally present at a given frequency of 



excitation. Above the plasma frequency, both the trans- 
verse and longitudinal modes are propagating, while be- 
low the plasma frequency both modes decay exponen- 
tially. In considering boundary value problems, it is clear 
that a wave incident on a half space filled with a nonlocal, 
plasmonic medium will generally couple to both types of 
waves. To avoid the system being underdetermined, an 
additional boundary condition must be used as will be 
discussed in the subsequent section. 
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FIG. 3. A simple interface between a dielectric and a metal. 

Let us now consider a multilayered structure, that 
could be as simple as the single interface shown in fig- 
ure [31 invariant in two directions, here taken as x and y. 
The z axis is thus perpendicular to any interface consid- 
ered, as shown in Fig. [3] Without any loss of generality, 
it is possible to assume solutions that are translation- 
ally invariant along the y (out-of-plane) direction. As 
shown in the appendix, the system of equations (jHJ and 
(|10[) can be split into two subsystems corresponding to 
s (electric field polarized perpendicular to the plane of 
incidence) and p (magnetic field polarized perpendicular 
to the plane of incidence) polarizations. Moreover, we 
will assume from now on that all the fields present an x 
dependence that varies as e lkm x (or, equivalently, we take 
the Fourier transform along the x axis). 

For the s polarization, V x E = yields E y = 0, so 
that no bulk plasmon can be excited. Nonlocality has 
then no impact on this polarization, so that we will deal 
in the following with p polarization only. 

Equation (|16l) then yields 

d 2 z H y = -{ekl-kl)H y , (22) 
so that the magnetic field can be written 

H y = (A e ik >* z + B e- ifc « z ) e ik * X ~' IUJ ' (23) 

with k zt = \/ekQ — k 2 where ko = — . The E x and E z 
accompanying fields can be found using equations 

E x = -^—d z H y (24) 

E z = --^—d x H y . (25) 
iuicq e 

For longitudinal waves, the wave equation (fH?)) be- 
comes 

8 2 Z E X = (kl + -) E x . (26) 
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For normal incidence (for k x — 0), depending whether 
u) is smaller or larger than . the bulk plasmon 

will be respectively evanescent (Sft(a) > 0) or propaga- 
tive (3?(a) < 0) . In the visible range, we usually have 
uj < ,f p so that the bulk plasmon is evanescent and 

the above equation can be solved to yield 

E x = {Ce K ' z +De- KlZ )e lk * x - lLUt (27) 

with 

The fact that longitudinal waves are curl-free yields 

E z = ~^d z E x , (30) 

which allows determination of the contribution of the 
bulk plasmon to E z if needed. 

IV. ADDITIONAL BOUNDARY CONDITIONS 

The nonlocal nature of the metal results in the appear- 
ance of a longitudinal bulk plasmon mode that can be ex- 
cited from the metal interface, in addition to the surface- 
localized plasmon polariton. The well-known Maxwell's 
boundary conditions are not sufficient to uniquely de- 
fine the amplitudes of these independent waves. More 
specifically, for each metallic layer, two new unknowns 
are introduced and must be resolved in the solution of the 
electromagnetic boundary value problem. To avoid deal- 
ing with an underdetermined problem then, additional 
boundary conditions must be imposed at the metal inter- 
face. 

The issue of boundary conditions has been abun- 
dantly discussed in the context of spatially dispersive 
crystals, and a variety of different boundary conditions 
has been propose d 44 ' 47 . In the context of the hydrody- 
namic model, the choice of boundary conditions is much 
simple r 20 ' 22 , essentially because fewer types of waves are 
involved. Two additional boundary conditions are typ- 
ically considered in the case of an interface between a 
metal and a dielectric, when the contribution of the 
bound electrons is overlooked: either (i) P z = 20 ' 22 or 
(ii) the continuity of E z 14 ' 48 . If the considered dielec- 
tric is vacuum, then these two conditions are equivalent. 
Condition (i) can be justified because the polarization in 
the metal is due to actual currents; since the electrons are 
not allowed to leave the metal, then the normal current 
must vanish at the interface and also the polarization. 
Condition (ii) can be justified by treating the interface 
as smooth for all fields, including the normal component 
of the electric field. 



In our description of the response of metals, the sus- 
ceptibility attributed to bound electrons, \b is consid- 
ered purely local. One might expect that the equation 
V • D = would impose a supplementary condition 
(namely the continuity of D z ) 1 leaving no freedom in the 
choice of the boundary condition. This is however not 
the case: in multilayered systems, the continuity of H y 
through an interface implies the continuity of D z , so that 
an additional boundary condition is still required. 

The response of the metal in our description is partly 
the response of a standard dielectric medium, so that 
there is no reason to assume the continuity of E z at 
the surface of the metal. Condition (ii) thus appears 
very difficult to support when the contribution of bound 
electrons is taken into account as a local, polarizable 
medium. 

The underlying physic o 20 i 22 behind condition (i), that 
free electrons cannot escape the metal, does however not 
lead here to P z = at the edge of the metal because 
not all the polarization comes from actual currents in 
the free electron gas. It is thus not reasonable to use 
boundary condition (i) for the case when bound electrons 
contribute to the polarization response. 

It would be physically reasonable to consider that only 
the polarization linked to actual current leaving the metal 
should be zero at the interface between a metal and a 
dielectric. For multilayered structures, this condition can 
be written 

Pf, = 0- (31) 

at the interface as an additional boundary condition. We 
underscore that this boundary condition is not equivalent 
to conditions (i) and (ii) when the outside medium is 
vacuum. 

In the case of an interface between two metals, again, 
condition (i) is hard to justify, but the interface obviously 
should not be considered as impervious to free electrons. 
Instead, it would sound to consider that the currents, 
and thus the polarization P/, should be continuous. This 
would actually provide the two additional boundary con- 
ditions that are required for an interface between two 
metals. Although we will not consider here structures in- 
volving such an interface, we emphasize that taking into 
account the contribution of bound electrons to the re- 
sponse of metals seem to lead to unambiguous boundary 
conditions based on physical reasoning. 



V. REFLECTION FROM A METALLIC 
SURFACE 

Let us now consider an incident plane wave coming 
from above (z > 0) and propagating in a dielectric 
medium with a permittivity e^, reflected by a metallic 
interface located at z = 0, as shown in Fig. [3]- the metal 
being characterized by a permittivity e. 

For p polarization, the magnetic field in the dielectric 



6 



region can be then written 



Hy = { 



— fv-ikz Z 



+ r e 



ik z z\ 



ik x x—iu) t 



(32) 



where k z = y/eak'o — A;^ and ko = — , while the electric 
field along the x direction has the form 



E x — 



ik z 



-ik z 2^ 



ik~ x—iuj t 



tuje e d 

In the metal, the magnetic field can be written 

jj z ^ik x x — icu t 



where Kt = — e k^ , and the electric field 

™t a z , j-> Ki z \ ik x x—iu)t 



E r 



E, 



-Ae K 



-Ae K 



Be KlZ 
hi 



Be K 



Ak x x—iuj t 



(33) 
(34) 

(35) 
(36) 



The magnetic field H y and the x component of the 
electric field E x are continuous at z = so that 



1 



= A 



(r-l) l ^L = ^iA- 



zwep B. 



(37) 
(38) 



Since Pt = —-^d x H y — e (l + Xb)E z , the condition 
Pf = in the metal at the interface, can be written 



tk^r A I — 

e 1 + 



Xb 



= —j — iuj£o B 



Finally A and B can be eliminated to yield 

L - n 



trf 



where 



trf 



9. 



1 



Xb 



(39) 



(40) 



(41) 



The reflection coefficient indicates that the bulk plas- 
mon is not excited at normal incidence for k x = be- 
cause in that case only one component of the electric 
field is present in the incident and reflected fields. When 
the angle of incidence increases the excitation of the bulk 
plasmon is more and more important because of the in- 
creasing E z component. Of course ki is increasing too, 
which means that the bulk plasmon penetration is more 
shallow, but only slightly - so that the increases essen- 
tially as k\. 



with k z = \/kx — Ed &q = — ik z , then it is meaningless to 
define a reflection coefficient (|40|) . but still we can write 
that 



D 

C 



«a _ St +Q, 

£rf g 

Ssl + St _ Q 



(43) 



The surface plasmon is a solution for which D ^ and 
C = 0, thus corresponding to a pole of the left hand side 
of equation (|43l) . and a zero of its denominator, so that 
the dispersion relation can be written 



e 



(44) 



The larger the propagation constant k x , the larger Q 
and thus the larger the impact of nonlocality. However, 
for surface plasmons, very large values of k x are diffi- 
cult to reach (typically, the maximum effective index is 
around 1.4 for a silver-air interface) as shown in figure 
31 The impact of nonlocality on bare surface plasmons is 
thus very small. In the following, we will see that for a 
metal-dielectric-metal waveguide with a very thin dielec- 
tric layer, the impact of nonlocality on the guided mode 
is much more important because very large k x values can 
be reached whatever the wavelength. 
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FIG. 4. Dispersion relation (uu as a function of k x ) for a sur- 
face plasmon at the interface between silver and air. The 
local description (thick solid curve) can almost not be dis- 
tinguished from the nonlocal description (thick dashed curve, 
almost identical with the thick solid curve). In order to illus- 
trate the effect of nonlocality, we show here (dotted line) the 
impact of an exaggerated nonlocality (/3 multiplied by ten). 
The thin solid curve is ui — k x c. 



VI. SURFACE PLASMON 

If the field is not propagative in the dielectric region, 
but has the form 



H y = (Ce K > z + De~ K * 



z z\ -ik x x — iuj t 



(42) 



VII. THE IMPACT OF BOUNDARY 
CONDITIONS 

The form of the reflection coefficient and to the disper- 
sion relation above clearly show that 51 is the parameter 
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controlling the influence of the nonlocality on propaga- 
tion phenomena. Moreover, it makes manifest the conse- 
quences of a change in the boundary conditions. 

In the literature, the entire metal response is often at- 
tributed to the free electrons, while the response of bound 
electrons is neglecte d 20 ' 22 ' 48 . When the bound electron 
response is neglected, the dispersion relation of the bulk 
plasmon yields 



ki = 




instead of 



Kl 




(45) 



(46) 



where a local contribution from bound electrons is 
assumed 1 ^. 

If the boundary condition with the dielectric is chosen 
to be the continuity of the component of the electric field 
normal to the interface, then we have 



<> : — I - - ~ 



1 

e<2 



(47) 



where Ki can be calculated using one of the above ex- 
pressions, depending on the description of the metal's 
properties. When the entire polarization P z is chosen to 
vanish at the interface, we have instead 



Kl 



1 



(48) 



In the following, we will investigate all the different 
descriptions that are presented in table U to show that, 
even if they differ regarding the quantitative impact of 
nonlocality, they all at least agree qualitatively. 



Descr. 



A.B.C. 



K + - 

kl + j 
k 2 + ^ 

2 

U 2 J P 



+ 1 



1 + X6 



Xf 

i 

Xf 1 + X6 
1 

1 + X6 



E z continuous 

P 2 =0 
E z continuous 

Pfz=0 



-1 



l _ i 

i l+Xb 



TABLE I. Summary of the different descriptions of nonlocal- 
ity. The first two do not consider separately the contribution 
from the bound electrons, the last three do. The last one is 
the one that is preferred in this work. 



VIII. METAL-DIELECTRIC-METAL 
WAVEGUIDE 

While the impact of nonlocality can be considered mi- 
nor for the single interface problem above, nonlocal ef- 
fects can be far more evident in multilayer systems. In 



metallo-dielectric layers, it is possible to reduce the thick- 
ness of layers to the nanometer or even sub-nanomctcr 
scale; modes that propagate in such layers can be sig- 
nificantly confined, to the point where local models are 
forced to break down. For this reason, multilayer systems 
and structures based on multilayers can be useful as an 
experimental tool to investigate and measure nonlocal 
effects. 

In this section, we consider the case of a dieletric with a 
permittivity sandwiched between two metallic surfaces 
(as shown in figure [S]) and study more thoroughly the 
influence of the nonlocality of the metal on the first even 
guided mode (the fundamental mode). 



h/2- 




Dielectric 




FIG. 5. Metallic waveguide of width h. 



A. Dispersion relation 

We consider here a symmetric waveguide, the metal 
being the same on both sides of the dielectric layer. The 
magnetic field in the dielectric can be written as 

\ Jik x x—iuj t 



By = (Ce K * z +De- 



(49) 



As we have seen in the previous section, at z = — |, 
we have 



D e Kz 2 



s*. - £1 + 
«t -I- St -Q 



(50) 



while for z = +-| (the z axis has to be reversed, which 
means C and D should be exchanged) 

Ce*-4 s^-st+n 



D, 



2*. + St - Q- 



Combining these two equations, we get 

/ 2t - St +0\ 2 

and finally either the mode is symetrical (C 
we have r = e Kz which can be written 



— tanh 

t d 2 



Hi 

e 



(51) 

(52) 
D) and 

(53) 
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or the mode is antisymetrical (C = —D), which means 



a K. z h 



and finally 



f^Z i ^ Z fa* 

— coth 

e d 2 



(54) 



B. Nature of the guided modes 

We first discuss the nature of the guided modes in a 
thin metallic waveguide. There are two situations that 
are clear and for which the guided modes of the structure 
have well-posed definitions^ : 

• The perfect metallic waveguide, which supports a 
fundamental mode that is flat and that has no cut- 
off (it is supported whatever the thickness of the 
metallic waveguide). In addition, we have analyt- 
ical expressions for the propagation constant and 
field profile of all the modes. For the fundamental 
mode, we have 



(55) 



The plasmonic (i.e. wide) metallic waveguide, 
which supports coupled surface plasmons^. At a 
given frequency and for a wide enough guide, the 
even and the odd surface plasmon modes present 
propagation constants that can be arbitrarily close 
to the propagation constant of the surface plasmon 



e<2 



e( W )' 



(56) 



even for complex values of e(w). 



For the case of a thin (a few nanometers) waveguide 
we seek the best description to retain for the only guided 
mode found. 

Consider the case of coupled surface plasmons first. We 
can approach the condition of a perfect metallic waveg- 
uide by making the permittivity of the metal change such 
that its real part tends towards infinity. As can be seen 
in figurc[B]the odd mode tends towards the fundamental 
mode but the field inside the dielectric (index of 1.58) al- 
ways stays evanescent. The even mode tends towards the 
first even mode of the perfect metallic waveguide and the 
field becomes propagative at some point (where the real 
part of the propagation constant becomes smaller than 
the optical index of the dielectric). The point at which 
the field of the even mode becomes propagative could 
even be defined as a limit between the "coupled surface 
plasmon" and the "perfect metallic waveguide" pictures. 

Now consider starting with a large waveguide (500 nm) 
and decreasing its width down to a few nanometers. As 
can be seen in figure[71 the even mode presents an increas- 
ing propagation constant. The odd mode, by contrast, 
presents a decreasing propagation constant - the field in 




1 1.5 

Re(k x /ko) 



FIG. 6. Trajectory in the complex plane of the quantity for 
three different waveguided modes (solid line: first even mode; 
dashed line: first odd mode ; dotted line: second even mode) 
when the permittivity of the metal goes from e = — 10 + i 
(circles) to e = —2000 + i, for a thickness of the waveguide 
of 500 nm and a dielectric with an 1.58 optical index. The 
intermediate value of —100 + i is indicated on the curves. 



the dielectric even becomes propagative as in the previ- 
ous case. For thin layers (smaller than 128 nm here) this 
mode presents a very large imaginary part and a very 
small real part: it can be considered as evanescent in the 
x direction even if the cut-off cannot be defined precisely. 
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FIG. 7. Trajectory in the complex plane of the quantity 
for three different waveguided modes (solide line: first even 
mode; dashed line: first odd mode ; dotted line: second even 
mode) when the width of the metal goes from 380 nm (cir- 
cles) to 10 nm. The permittivity of the metal is taken equal 
to the permittivity of gold at 608 nm, e ~ —10.01 + 1.44i. 
Intermediate thickness of 221 nm (when the real part of the 
index of the odd mode becomes smaller than 1.58, so that 
the coupled plasmon picture becomes less relevant), 128 nm 
(when the odd mode can be considered non-propagative, and 
the coupled plasmon picture is not relevant any more) and 50 
nm are indicated on the curves. . 

For a small dielectric thickness, the waveguide thus 
behaves much more like a perfect metallic waveguide (a 
fundamental mode with no cut-off, no propagative even 
mode) and except for the fact that the field of the first 
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even mode is evanescent in the z direction, has not much 
to do with the coupled surface plasmons situation. This 
is why we refer to this mode as the fundamental mode of 
the waveguide. 

This mode is however sometimes called gap-plasmon in 
the literature^ 5 -, a term that underscores the differences 
between the actual mode and the fundamental more of a 
perfect metallic waveguide. 



C. Nonlocal effects 

When the waveguide becomes extremely thin, as can 
be seen in figure [7J the effective index (and thus k x ) of 
the fundamental mode (with a dispersion relation given 
by (|53p) can become arbitrary large. When k x is larger, 
n is larger, which means that the non-locality has a much 
larger impact on the mode's propagation constant. It is 
possible to compare (see figures 13 [9] for a waveguide filled 
with a dielectric with a 1.58 optical index) the local ef- 
fective index as a function of the waveguide's width with 
a local and with a nonlocal theory. Obviously the impact 
of nonlocality is limited for h > 5 nm but it can become 
very important under that threshold. The parameters we 
have considered for gold are given in^ and (3 = 1.27.10 6 
m/sS 



1000 
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FIG. 9. Same as figure [8] except for the log scale. 



the medium is nonlocal or not, as described in section 
IVIII and summarize in table Jj. 

The hydrodynamic model is often said to exaggerate 
nonlocal effects. It could be expected that taking into 
account the response of the bound electrons, which can 
be considered as local, would lower the impact of nonlo- 
cality on the guided mode compared to when the whole 
response of the metal is considered nonlocal. Figures [8] 
and [5] show that this is paradoxically not the case when 
the boundary condition that we consider as being the 
most physical (Pf = 0) is not chosen. Considering a 
separate response of the bound electrons actually lowers 
the effective plasma frequency, as explained above, which 
leads to a deeper penetration of the field corresponding 
to the bulk plasmons, which may in turn increase the 
importance of this held (depending, of course, on the 
boundary conditions). 

When the condition we propose is used, the impact of 
nonlocality is even lower than when considering a com- 
pletely nonlocal response of the metal and using P z = 0. 
This actually makes us think the boundary condition we 
propose here, is not only the most sound physically, but it 
may even yield an more accurate estimate of the nonlocal 
effects. 



FIG. 8. Effective index of the guided mode at 600 nm, as a 
function of the dielectric width, h. The dispersion relations 
are shown for different descriptions (see table HJ: the com- 
pletely local case (thin black curve, top), -P z (0) = (descrip- 
tion 1, dash-dotted line) and E z continuous (description 2, 
dashed line) with no identified contribution of the bound elec- 
trons, and descriptions separating the contributions of bound 
and free electrons, with P z — (description 3, dash-double 
dotted line), a continuous E z (description 4, thin dash line) 
and finally our description (description 5, thick solid curve), 
that is preferred in this work. 

Since different descriptions of nonlocality exist in 
the literature, we have compared our approach to the 
other descriptions available (different boundary condi- 
tions, as well as considering that the whole response of 



IX. CAVITY RESONANCES FOR METALLIC 
STRIPS COUPLED TO A METALLIC FILM 

Many structures and phenomena rely on the funda- 
mental mode of the metallic waveguide like the enhanced 
transmission by subwavelength slit arrays 4 ^—, highly 
absorbent gratings 5 - or strip nanoantennas^r— to men- 
tion a few. The latter are patches that are invariant 
perpendicularly to the plane (see Fig. [TU|). The mode 
that is guided between a strip and the metallic film, 
whose dispersion relation is given by (|53[) as long as the 
patch is thick enough, if reflected by the edges of the 
strip. The reflection coefficient r of the mode can be 
computed easily 5 ^ using a Fourier Modal Metho d 53 ' 54 . 
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When the strip is wide enough, Fabry-Perot resonances 
may occur—. 



200 nm 

d=50 nm^ 



h=3 nm I 




FIG. 10. Strips (rods with a 50 nm by 50 nm section) are 
separated from the metallic film by a 3 nm thick dielectric 
with an optical index of 1.58. The structure considered here 
is periodic, with a 200 nm period. 
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The local energy density (and hence the absorption) 
should be proportional to the square of the field ampli- 
tude, given by a Fabry-Perot formula—, yielding 



This model allows the accurate prediction of the posi- 
tion of the resonance when a purely local response of the 
metal is assumed, as shown figure ITT] the resonance pre- 
dicted using model (|57l) (k x being computed using dis- 
persion relation Eq. (|53l) with fi = 0) occurs exactly 
where there is a dip in the reflectance of the nanorods 
covered surface. This confirms the physical analysis of 
the structure and that a one-mode model is sufficient to 
describe the resonances. 

As we have shown above, when nonlocality is taken 
into account, the propagation constant k x of the guided 
mode differs from the purely local case. That is why the 
resonances of the nanorods can be expected to be very 
sensitive to nonlocality when the thickness of the spacer 
is typically smaller than 5 nm. 

Full COMSOL simulations based on the hydrodynam- 
ical model with the boundary conditions we suggest in 
this work (description 5 in table H| show that the res- 
onance of the structure is largely blueshifted compared 
to the purely local simulations. This is completely ac- 
counted for by model (I5T[) when using a propagation con- 
stant k x computed using the dispersion relation (|53]l and 
keeping the same coefficient reflection r as for the lo- 
cal case. This proves that nonlocality intervenes almost 
only through the change of the propagation constant of 
the guided mode, and not at all through a change of the 
reflection coefficient. 

Such structures, or structures presenting a very simi- 
lar behaviour—, are obviously a way to assess experimen- 
tally the effects of nonlocality on the guided mode of the 
metallic waveguide with a good accuracy. 

X. CONCLUSION 

We have proposed in this work an improvement of the 
hydrodynamic model by clearly separating the nonlocal 



FIG. 11. Bottom : Reflection spectrum according to local 
(RCWA, dashed line) and nonlocal (COMSOL, solid line) sim- 
ulations. Top : Model (|57[1 for the field intensity under the 
strips for the local (dashed line) and nonlocal (solid line) the- 
ory. The agreement between the simulations and the model 
are in excellent agreement for the local as well as for the non- 
local theory. 



response of the free electrons, and the response of the 
bound electrons, considered as locals. Such a distinc- 
tion makes the discussion about the additional boundary 
conditions much more clear, leaving nothing but a sin- 
gle condition that seems physically sound: no current of 
free electrons leaving the metal. We have shown that 
this condition leads to a lower impact of the nonlocal- 
ity than many other descriptions based on the hydro- 
dynamic model. This description may thus answer two 
main concerns regarding this kind of models compared 
to Feibelman's approac h 26 ' 28 : the uncertainty about the 
boundary conditions and a tendency to exagerate the ef- 
fects of nonlocality. Furthermore, recent experimental 
results have shown that the hydrodynamical model can 
describe nonlocal effects very accuratelyi 6 .. Given the re- 
duced complexity of the hydrodynamic model relative to 
full quantum and other microscopic models of electron 
response, it is of continued interest to further explore 
the accuracy of these models in the context of plasmonic 
nanostructures. 

Following previous work on slot waveguides that sup- 
port gap plasmons^i, we have shown that the slow light 
regime reached when the waveguide is only a few nanome- 
ters thick is responsible for a large enhancement of the 
nonlocal effects. Using these results, we have studied the 
impact of nonlocality on patch nanoantennas and shown 
that it should be easy to detect, paving the way for future 
experiments. 

Our analysis is of relevance to numerous nanophotonic 
devices, including metallodielectric waveguides, nanoan- 
tennas and nanocavities, which rely on the excitation of 
gap plasmons on very thin, conducting layers for their 
operation^— . These resonant structures have a variety 
of diverse applications, for instance, as highly efficient 
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concentrators and absorbers of light 17,52 . The descrip- 
tion of conductors we provide can also prove useful when 
testing the limits of the classical theory for describing 
structures containing metals or doped semiconductors, 
for which the response of the bound electrons are strong. 

As has been once more shown here, the hydrodynamic 
model yields analytical results that help to understand 
the underlying physics of nonlocalit y 24 ' 48 . It presents the 
supplementary avantage of being easy to use in simula- 
tions with complex geometries 16 . The analytical calcu- 
lations we have presented, beyond the clarification they 
may brin g 41 ' 42 , are thus a first step towards the exten- 
sion of widely used numerical methods^— to account 
accurately for nonlocality. 



written 



d z E x — d x E z — iuj /io H y 
1 



-ad£E x +E X - ad x 8 z E z = 



-ad z z E z +E Z - ad x d z E x = 



1 



lujeo e 



d z H y 



d x H y 



(64) 
(65) 

(66) 



By applying the operator —ad x d z to equation (f6"6l) . 
operator 1 — adl to equation (|65| and subtracting one 
resulting equation from the other, one gets 



(67) 



Repeating the same procedure, but applying 1 — ad x 
to equation (ITjBI . —ad x d z to equation (|65p . and subtract- 
ing the resulting equation from the other we obtain the 
decoupled system of equations 



APPENDIX I 



Let us write equations © and (ITTJ|) within the metal 
in Cartesian coordinates in the case where the fields do 
not depend on y : 



d z E x - d x E z = iuj /i H y (68) 
(1 - a {dl + dl)) E x = ^~ a ^ + 9 ^ d z H v (69) 

(1 - a (dl + dl)) E z = -^—^±31 dxHy , (70 ) 

We can apply the inverse of the differential operator 
1 — a (9 2 -I- dl) to both sides of the equations to obtain 
the classical system 



-d z E y = iui H x (58) 

9*^ - d x E z = iuj /i H y (59) 

d x E y = iuj fi H z (60) 

-d z H y = -iuj€ e (E x - ad 2 x E x - ad x d z E z ) (61) 

d z H x - d x H z = -iu!E eEy (62) 

d x H y = -iuje e (E x - ad 2 z E z - ad x d z E x ) (63) 



ad 2 z E z - : 



d z E x - d x E z = iuj /io H y 
1 



E x — — 

iujeo e 

E, = --± 



d z H y 



d x H y . 



(71) 
(72) 

(73) 



This system is identical with that corresponding to a 
purely local response of the metal and its solution satis- 
fies V • E = 0. The wave that it describes is referred to as 
the transverse wave because when it is propagative, the 
electric field is orthogonal to the propagation vector. 

But to this solution should be added any solution for 
which 



This system of equations can be split into two subsys- 
tems corresponding to s (electric field polarized perpen- 
dicular to the plane of incidence) and p (magnetic field 
polarized perpendicular to the plane of incidence) polar- 
izations. The s subsystem is identical to the subsystem 
without taking nonlocality into account, because of the 
simple form of equation (f6"2"j). Nonlocality has then no 
impact on this polarization, so that we will deal in the 
following with p polarization only. 

The subsystem concerning the p polarization can be 



-a (d 2 x + dl) E x + E x =0 
-a (dl + d 2 z ) E z +E z =0 
-a(dl + dl)H y + H y = 0. 



(74) 
(75) 
(76) 



because it would also be a solution of system [B5] Using 
(|74p along with (1551) and (|6~4"|) . it is not difficult to show 



that H y = so that this solution satisfies 



d 7 E x = d x E, 



(77) 
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